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Abstract: The inclusion of B-factory data in fragmentation function analyses will improve 
the precision with which these are known. To remove the large 6-quark contribution from on- 
resonance data, a cut on the thrust event shape T is used, leading to a dijet final state. Here 
we complete our analysis of unpolarized fragmentation of (light) quarks and gluons to a light 
hadron h with energy fraction z in e'*"e~ — )• dijet + /i at the center-of-mass energy Q = 10.52 
GeV. In addition to the next-to-next-to-leading order resummation of logarithms of 1 — T, 
we include the next-to-leading order (NLO) nonsingular 0(1 — T) contribution to the cross 
section, the resummation of threshold logarithms of 1 — z, and the leading nonperturbative 
contribution to the soft function. Our results can be directly used to study correlations 
between fragmentation and the thrust cut imposed on B-factory data, which must be taken 
into account to extract fragmentation functions with reliable uncertainties. These correlations 
are also observed in Pythia, but are surprisingly smaller at NLO. 
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1 Introduction 

Hard QCD processes give rise to highly virtual partons that evolve by emitting radiation, 
until they hadronize into the nonperturbative states we observe. The physical basis of this 
"fragmentation" has not been constrained as much as jet production or deep-inelastic scat- 
tering (DIS), making this an active area of research. A quantitative understanding of parton 
fragmentation clearly extends the set of reactions that can be handled within a perturbative 
QCD approach. Applications include, for example, the investigation of the spin structure of 
the nucleon in semi-inclusive DIS [1, 2] and the study of hadron production at high pj- in 
collisions which is crucial to determine more accurately the relative suppression of hadron 
spectra (jet quenching) seen in heavy-ion collisions [3]. Interestingly, for many of these ap- 
plications, additional cuts on the hadronic final state are required to reduce the background 
from other processes, to help identify underlying partonic structures or to have better control 
on the initial hard-scattering kinematics [1-6]. 

A key observation is that the fragmentation process involves physics at well-separated 
energy scales. It has been shown that high-energy processes with an observed hadron in the 
final state can be described by factorizing the short-distance (partonic) physics associated with 
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the hard scale Q, which is perturbatively calculable in QCD, from (universal) nonperturbative 
long-distance contributions, see e.g. ref. [7]. At leading power in Aqcd/Qj the (unpolarized) 
fragmentation functions (FFs) D^{z,n) [8, 9] encode the information on the nonperturbative 
transition from an energetic parton i = {g, n, u, d, . . . } to a hadron /i, which carries a fraction z 
of its energy^, plus a remainder X. The knowledge of both perturbative and nonperturbative 
ingredients in factorization theorems is crucial to obtain theoretical predictions. Unpolarized 
FFs, for example, serve as an input to extract the flavor separated helicity distributions from 
spin asymmetries observed in polarized semi-inclusive DIS [1-3]. 

Parameterizations for the FFs have been constrained by fitting to data for single-inclusive 
charged hadron production in e"*"e~ at next-to-leading order (NLO) in perturbation the- 
ory [10-12]. More recently, global analyses have been performed to incorporate also semi- 
inclusive deep- inelastic scattering and/or pp^ pp data from HERA, RHIC and the Teva- 
tron [13-15]. To illustrate the current level of accuracy, the dominant D^^{z,fi = mz) is 
determined with uncertainties at the 10% level for z > 0.5 [12]. The FFs of the gluon and 
the non-valence quarks are known even less accurately. 

The Belle collaboration is currently studying light-quark fragmentation [16]. Since most 
of the current e~^e~ data is at Q = mz, this will significantly improve our knowledge of 
FFs. In particular, it will give access to the gluon FF via mixing with the quark FFs in 
the renormalization group evolution. To remove the large 6-quark background in their high- 
statistics on-resonance data, the Belle collaboration imposes a cut on thrust, which is an 
event shape variable defined as [17] 

T = maxi^^, (1.1) 

where the sum is over all final-state particles. More precisely, it is found that a thrust cut 
ofr = l — T<0.2 removes 98% of the B data leaving the thrust distribution dominated by 
the fragmentation of light (n, d, s) and charmed quark pairs [18]. With this cut the hadronic 
final state is given by back-to-back jets. 

In refs. [19-22] we developed the theoretical framework to study the fragmentation of 
a light parton into a light hadron h inside an identified jet. With "identified" we refer to 
jets determined either through an event shape like thrust [19, 20] or via a jet algorithm [22]. 
Our analysis provides the theoretical tools for a more exclusive study of fragmentation than 
is currently done, which is suitable for processes where additional cuts on the hadronic final 
states are imposed, such as the previous Belle example [18]. In particular, our framework 
allowed us to calculate the non-trivial correlations between and the energy fraction z in 
the observed cross section. Due to these correlations, the thrust cut modifies the shape in z. 
As already shown in the preliminary study of ref. [20], where we resummed the logarithms 
a" In™ T to next-to-next-to-leading logarithmic accuracy (NNLL), the effect is expected to be 
sizable and important for the analysis of fragmentation in B-factory data. Here we improve 

'^Due to perturbative radiation, the z in the fragmentation function is not the same as the experimentally 
measured momentum fraction [see e.g. eq. (2.1)]. 
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Figure 1. Correlations between the thrust cut t <t'^ and the observed momentum fraction z in the 
cross section of e+e" — > dijet + tt"*" for Q = 10.52 GeV. Curves and bands are plotted relative to the 
case — 0.2. A cut on thrust changes the shape in z. Left panel: our NNLL+NLO results with 
perturbative uncertainty bands. The result at NLO (dotted lines) contains negligible correlations. As 
input we use the HKNS fragmentation functions Z?f [12] at NLO. Right panel: the same observable 
using Pythia. 

our study in ref. [20] by including the NLO nonsingular contribution to the cross section, the 
resummation of threshold logarithms oi 1 — z and the leading nonperturbative correction to 
the thrust soft function, which are discussed in sec. 2. 

The correlations we obtain from our new analysis are shown in the left panel of fig. 1 for 
various thrust cuts t^, including the theoretical (perturbative) uncertainties. As can be seen 
there, the effect of the thrust cut is much stronger in the small z region than in the large 
z region. This has a straightforward physical explanation: for larger values of z, most of the 
jet energy is carried by the observed hadron, causing the jet to be more collimated and less 
affected by the thrust cut. In particular, for z > 1 — t^, the cross section is unaffected by 
the thrust cut. On comparison with the right panel, we see that the effect we find is fully 
compatible with the outcome of the Pythia event generator based on LL parton showering 
and the Lund string fragmentation model ^. We stress that the plot in the left panel has 
been obtained from our purely perturbative calculation, taking only the phenomenological 
parameterization of the universal -Df^ in ref. [12] as input. There we also show the correlations 
that we obtain at NLO (dotted curves), which are much smaller. Presumably, the single 
additional emission at NLO is insufficient to reliably describe this doubly-differential cross 
section, compared to the multiple emissions in our resummed calculation and in Pythia. 
The purpose of our study is to pave the way for a reliable extraction of the light-quark 
fragmentation functions using high-statistics on-resonance B-factory data, with controllable 
theoretical uncertainties. 

The paper is organized as follows. In sec. 2 we briefly review our theoretical framework, 
where we describe the updates of our analysis compared to our previous work in ref. [20] . In 

■^The results for other light hadrons, like the proton, are very similar to the ones shown here for tt^. 
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Figure 2. A schematic of the various subprocesses in e+e" — dijet + h: The (green) vertex denotes 
the hard interaction H, the (blue) jets are described by J and J'ij (the latter for the jet in which the 
hadron is observed), the fragmentation j — > /i is described by the fragmentation function and the 
effects of the (orange) soft radiation are contained in S. 

sec. 3 we discuss the features of the z-spectrum and r-spectrum and illustrate the effect of the 
different ingredients of our calculation. In sec. 4 we propose a first, straightforward, quanti- 
tative analysis of correlations between z and t'^ in Belle data. All the necessary theoretical 
input for this moment-space analysis is collected in app. C. App. A contains the details of the 
calculation of the NLO nonsingular contribution to the cross section. A detailed discussion 
of our choice of renormalization scales, as well as the scale variations used to estimate the 
perturbative uncertainties, is given in app. B. 

2 Calculating correlations between thrust and z in e+e — dijet + h 

We focus on the case of spin-averaged fragmentation in e~^e~ — t- dijet + h where one restricts 
to the dijet limit by a cut on thrust. Starting from the well-known factorization theorem 
for the inclusive measurement of thrust [23-26], in refs. [19, 20] we obtained a factorization 
theorem for the observed cross section at leading power in Aqqb/Q, which has this form: 

E /'^ c,(.si,cV)^ (.1, 

j=g,u,u,d... 

Q is the center-of-mass energy. The Dj{x) denote the standard, unpolarized fragmentation 
functions, which describe the hadronization j — )• h(x) at leading power. The index j runs 
over all parton flavors, including the gluon. The coefficients Cj are calculable in perturbation 
theory and contain double logarithms of t^, a" In'^r'^ (m < 2n), and "threshold" logarithms 
of 1 — z. Our setup enables us to resum all these large logarithms in the dijet region of the 
thrust distribution, which is relevant for B-factories. This is necessary to achieve perturbative 
convergence with small uncertainties in Cj when < 0.2 and when z > 0.5, as we will show 
in sec. 3. The scale fj, in eq. (2.1) is arbitrary, but must be chosen to be the same in Cj and 

We obtained the expression for Cj by exploiting factorization in the framework of soft- 
collinear effective theory [27-30]. The underlying idea is that dynamics at well-separated 
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scales combine incoherently, allowing one to factorize their contributions to the cross section. 
In the dijet region, the relevant scales and the subprocesses they correspond to are illustrated 
in fig. 2 and eq. (2.2). Specifically, we distinguish the scale Q for the hard collision e+e" — )• qq, 
the scale ^/t^Q for the jet production (showering), the scale t^Q associated with the soft 
radiation between the jets, and the scale Aqcd at which hadronization takes place. The 
coefficients Cj(T^,Q'^,x/z, fi) are given by the convolution of functions encoding the effects 
at these different perturbative scales plus a nonsingular piece which contains contributions 
that do not factor in this way. Schematically, 

X 

z ' 



q 



(2.2) 



where cjg is the tree-level cross section. Since it is not known whether the observed hadron 
h fragmented from the quark or the antiquark jet, we sum over both possibilities. The 
C"^ contains the terms of the NLO cross section that are not enhanced by logarithms of 
r, and is 0{t) suppressed. The calculation of the nonsingular contribution to d'^a^ /{dr dz) 
at NLO is discussed in detail in app. A and its numerical effects are shown in sec. 3. All 
the other ingredients in the Cj's have been calculated in the literature to some order in the 
perturbative expansion. We refer to [20] for details. In particular, our one-loop determination 
of the coefficients J7ij in ref. [20] (partly contained also in ref. [31]) enabled us to analyze the 
cross section in eq. (2.1) up to NNLL accuracy. We stress that since fragmentation takes 
place inside a (hemisphere) jet of invariant mass of order t^Q"^, the momentum fraction z 
cannot be too small, to avoid contributions from soft hadrons outside the jet. The hadron 
momentum must be larger than the soft scale, implying z > t^. 

The factorization in eq. (2.2) enables the resummation of logarithms of r and 1 — z, as 
well as large tt^ terms in H, which is accomplished by evaluating H, J'ij, J and S at their 
natural scales and using their respective RGEs to evolve them to a common scale fi. We 
have made improvements in our choice of scales compared to our earlier work in ref. [20], 
specifically, we resum threshold logarithms and take the upperbound t < 1 — z into account. 
A detailed discussion of our choice running scales is deferred to app. B. 

We now comment on the various corrections in eq. (2.2). There are corrections of order 
Aqqj-,/[(1 — z)t^Q'^] associated with the ^Tij's. These coefficients describe the emissions from 
the parent parton (here a quark) at large virtualities building up the jet in which the hadron 
is identified. The associated corrections get sizable if the jet invariant mass gets small. In our 
plots in sec. 3 we have included the leading nonperturbative contribution to Sr, according 
to the analysis in ref. [32], which pushes the nonperturbative corrections to the soft function 
to order AQQj-,/(r'^(5)^, as indicated in eq. (2.2). These nonperturbative corrections become 
large for small t^, up to about 15% for = 0.15. We therefore need to be not too small. 
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Quark and hadron masses are treated as negligible in our calculation, and we briefly 
comment on the size of the corresponding kinematic corrections. The corrections due to the 
hadron mass m/^ are of order mf^/ {z^Q'^) [33]. For light hadrons, like pions and kaons, they 
are indeed negligible in the region where our framework can be applied. In the plots of sec. 3, 
we restrict q in eq. (2.2) to be tt, cZ, s. The correction due to the charm mass is expected 
to be of order m^/ {t'^Q'^) when charm is a valence quark, which is less significant than the 
other corrections to eq. (2.2). If the hadron h does not contain a charm valence quark, the 
only contribution to eq. (2.1) that does not get further suppressed by the smallness of the 
corresponding FF, is the one involving D^. 

3 Numerical analysis for e^e~ — dijet + vr"*" 

Here we illustrate our key results for the cross-section in e~'~e~ — )• dijet -l-vr^, given in eq. (2.1). 
As input we use the HKNS fragmentation functions D^^ and D'^^ [12] at NLO. To be consis- 
tent with ref. [12], we set as{^J. = rnz) = 0.125 with two- loop running, matching continuously 
across the quark thresholds In our analysis we do not include the effects of the uncertainties 
associated with as{mz) or with the FF parameters. Therefore we stress that the bands shown 
in our plots correspond only to the perturbative uncertainties obtained from scale variations, 
as explained in app. B. 

We have made the following improvements with respect to the plots shown in ref. [20]: 

(a) We include the NLO nonsingular contribution, which we calculate in app. A. This small 
0{t) correction plays a role for r"^ > 0.2. 

(b) We resum the threshold logarithms oi 1 — z according to our ref. [22], where we intro- 
duced a joint resummation of these logarithms and those in (suitable for values of 
in the range of interest for Belle). As described in app. B, this is accomplished through 
the choice of the renormalization scale for Jij. These effects are important for z > 0.5. 

(c) We include the leading nonperturbative correction to Sr according to ref. [32]. This 
amounts to a shift in the thrust distribution, r'^ — )• r'^ — 2Vii/Q, where Cli = 0.252 GeV 
in the MS scheme. The effect of this is most prominent at small t'^. 

In figs. 3 and 4 the effects of the thrust cut and our resummations are illustrated for fixed 
values of z. As the cut becomes stronger (i.e. for smaller t^) the cross section gets reduced, 
and by a larger amount in the case of our resummed calculation. For < 0.2 our resummed 
result starts to differ from the leading order (LO) and NLO results, showing the effect of the 
resummation of double logarithms of t'^ and the leading nonperturbative correction. 

The nonsingular NLO correction is singled out in fig. 4, where we switch off both the 
threshold resummation and the resummation of vr^-terms in the hard factor H of eq. (2.2). 
The nonsingular contribution is not large, but it causes our resummed calculation to merge 

^In the next sections, which do not involve HKNS FFs, we wiU use a more reahstic value for as- 
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Figure 3. The effect of the thrust cut in e+e^ dijet + 7r+ for Q = 10.52 GcV and z = 0.5 at LO, 
NLO and NNLL+NLO. The curves are normalized at t'^ = 0.5. The bands show the perturbative 
uncertainties obtained from scale variations discussed in app. B. 




T T 

Figure 4. The effect of the inclusion of the NLO nonsingular contribution is illustrated for z = 0.5 
(left panel) and z = 0.8 (right panel), with tt^ and threshold resummations switched off. 



with the NLO in the r'^ — )• 0.5 region where resummation is unimportant, and removes 
an unphysical decrease of the integrated cross-section for large values of r^. If we include 
threshold and vr^ resummation, we should also include these resummation in the nonsingular 
contribution to the cross section, to preserve the above properties. Following ref. [34], we 
estimate the vr^ resummation of the nonsingular (which is formally higher order in as) by 
simply multiplying it with the corresponding evolution of the hard function UniQ"^ , ^j-h, I/^hI), 
given in eq. (C.7) of ref. [20]. We have not attempted to resum the threshold logarithms in 
the nonsingular, so our calculation does not (yet) reliably describe the region where both r 
and z are large. 

Next, we look at the effect of the thrust cut on the z spectrum ("the shape of the 
fragmentation function"), which is shown in figs. 5 and 6 for the default value t'^ = 0.2. 
In fig. 5 we compare our resummed result at next-to-leading logarithmic order (NLL) and 
NNLL+NLO with and without threshold resummation. The threshold resummation leads 
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Figure 5. Comparison of NLL and NNLL+NLO with/without threshold resummation, for r'^ — 0.2. 




Figure 6. The effect of the thrust cut t'^ = 0.2 on the z-spectrum for e+e^ — > dijct + 7r+ with 
Q — fO.52 GeV at LO, NLO, NLL and NNLL. The right panel shows the curves and bands in the left 
panel relative to the central curve for the LO result. 



to slightly better convergence and smaller uncertainties at NNLL+NLO for z > 0.5. The 
effect is not as dramatic as in ref. [22] because the scales are much smaller here. In fig. 6 we 
compare our NLO, NLL and NNLL+NLO curves with the LO result, which is independent 
of T^. The differences between the shape in z are not entirely due to the thrust cut, but also 
correspond to different orders in (resummed) perturbation theory. As is clear from fig. 1, we 
also find correlations between and z. These effects have to be taken into account when a 
fit to experimental data is performed to extract the FFs. 



4 Studying t'^-z correlations in Belle data 

In fig. 1 and sec. 3 we showed that the thrust cut has a sizable, calculable effect on the spec- 
trum of the fragmentation variable z. Our framework allows one to study this qualitatively 
and quantitatively in the Belle data. 
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For a first quantitative analysis, we propose working in moment space, where eq. (2.1) 
simplifies to 

/■I A^h 

dzz^-i— (r^) = ^^,^(iV,^)C',(r^iV,Q,/.), (4.1) 

with 

M{N,n)= [\xx^-^D^{x,f^). (4.2) 
Jo 

At variance with eq. (2.1), we use the notation j = q,g by defining 

^' = Q«(^' + ^^ + ^c+^5) + Qd(^d + ^j + ^' + ^s-), c, = ajQl^Ca/Ql, (4.3) 

with Qu = 2/3 and Qd = —1/3 the electric charges of the up- and down- type quarks. In 
app. C we collect our numerical values for Cj(r^, N,Q = 10.52 GeV, = 1 GeV) in nanobarns 
for j = q,g, t"" = 0.15, 0.16,..., 0.2 and = 3, 4, . . . , 8, at NNLL+NLO. We use a,(mz) = 
0.118 and three-loop running ^, and employ two-loop splitting functions [35-37] to evolve 
to /i = 1 GeV. We provide a different set of coefficients for each of the independent scale 
variations in app. B, which allows one to take correlations between the uncertainties into 
account. 

Apart from the different electric charges for the up- and down- type quarks, the difference 
between Cu and is less than one percent. Thus the thrust cut does not provide a way to 
extract FFs for different quark flavors, but gives access only to the combination in eq. (4.3). 
Varying the thrust cut allows one in principle to separate out the gluon contribution. How- 
ever, at the scale = 1 GeV, where the fragmentation functions are usually extracted, this 
contribution to the cross section is small and dominated by uncertainties. 

Taking as input the Belle spectra, where both r'^ and z are varied, the numbers in tables 
2 and 3 can be directly used to perform fits for the moments Dq{N,fi = IGeV). If these 
fits produce acceptable values of x^/d.o.f., the next step will be the more involved analysis 
of correcting the measured z-spectra for the effect of the thrust cut, to reliably extract the 
FFs. Note that Mellin moment methods have been used in analyzing fragmentation functions 
before [13, 14]. We stress that our formalism enables one to use data at different r'^ to better 
constrain the parameters for the FFs. 

Experimentally, using moments has the disadvantage that they can be sensitive to the 
problematic z < 0.2 region, and the z > 0.9 region, where the experimental uncertainties 
are large [16]. The relative contribution of these regions to the moments of D^^ and is 
shown in fig. 7. This is why we restricted ourselves here to moments in the range 3 < A < 8. 
In addition, the thrust cut should neither be too soft (contamination with b quarks) nor too 
strong (reduced signal). As we explained in sec. 2, the requirements for our calculation to be 
under theoretical control are similar. 

*To avoid crossing flavor thresholds in our calculation, we choose to move the charm and bottom thresholds 
outside the range of scales that we work in. 
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Figure 7. The composition of the Nth moment of (left panel) and (right panel). Using 
different colors we show the contribution from the problematic regions z < 0.2 and z > 0.9. The 
contribution from z > 0.9 for d quarks is very small and hardly visible. 

5 Conclusions 

We have studied the cross section for fragmentation with a thrust cut, which is relevant to 
analyze fragmentation in high-statistics on-resonance B-factory data. The shape of the cross 
section in z is altered beyond LO and includes a r-dependent contribution. The correlation 
between r and z found in our resummed calculation is observed in Pythia as well, but is 
much smaller at NLO. Our plots in sees. 1 and 3 and the discussion in sec. 4 provide the 
starting point for a qualitative and quantitative analysis of this effect. We emphasize that 
binning the data in allows one to study the correlations in and z, and leads to more 
stringent fits for the fragmentation functions. 
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A Nonsingular contribution at NLO 

We obtain the NLO nonsingular correction to the cross section for e~^e~ — >• dijet + h, by 
subtracting the singular In™ t'^ terms from the full NLO calculation. The singular NLO can 
directly be obtained from our resummed calculation in ref. [20], so we start by calculating 
the full NLO cross section. Although most ingredients are well-known, we include them here 
for completeness. 
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A.l One-loop diagrams and kinematics 

We perform the calculation using dimensional regularization for the UV and IR divergences 
in the MS scheme. The tree-level cross section e'^e~ — )• qq is given by 



(0) 



(A.l) 



The real and virtual contribution at one loop are [38] 
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We employ the usual parametrization for the e^e 

'ipi-q CMS '^Pi 



Xi 



> qqg final state 
< < 1 , xi + 3:2 + X3 = 2 , 



(A.3) 



where q is the momentum of the e+e~ pair and the Pi=i,2,3 are the momenta of the quark, anti- 
quark and gluon, respectively. The phase-space integration boundaries for the real emission 
(jj^^ in eq. (A.2) depend on the observable. These are quite simple for the observables we are 
interested in: the Xi are the momentum fractions and thrust T = 1 — r = maxjxi, X2, 2:3}. 

A.2 Quark and gluon fragmentation 

We first consider the cross section where only the momentum fraction z of the hadron h is 
observed. This is described by the factorization theorem 
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up to Aqcd/Q power corrections. Replacing /i by a quark or a gluon, we find at one-loop 
order that 
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restricting to one quark flavor for simplicity. Here we used that 



da, 



(0) 



do", 



(0) 



0, 



da; dx 
as well as the partonic fragmentation functions in pure dimensional regularization, 



(A.6) 



= 5(l-z) + 0(a,). 



(A.7) 



The one-loop splitting functions are [39] 
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Measuring the momentum fraction of the (anti)quark or gluon, the real radiation contribution 
in eq. (A. 2) becomes 
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The virtual piece in eq. (A. 2) gives a contribution ay^6{l — z) to (t^\ which largely cancels 
the divergences of Cg^jj- The remaining infrared divergences are exactly those required by 
eq. (A. 5), and the finite terms yield 
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+ 



Inz 
1 - z 



3 1 



4(l-z) 



+ 



da. 



(1) 



dz 



a 



vr^ 9 
3" ~ 4 

(0) OLsCp 



vr 



Pgg(z) -ln^+ln(l-z) + 21nz 



(A.IO) 



This agrees with the calculations performed a long time ago in refs. [36, 37, 40]. 
A. 3 Fragmentation with a thrust cut 

The easiest way to calculate the fragmentation cross section with a thrust cut at NLO is 




Figure 8. The cross section for e+e^ dijet + tt"*" at Q = 10.52 GeV, differential in the momentum 
fraction z and thrust t, and separated into its singular and nonsingular contribution. Since these 
separate pieces can be negative, absolute values are plotted. 



The first term on the right-hand side was calculated in the previous section. The second term 
is described by a factorization theorem very similar to eq. (A. 4) 



do-fc 
dz 



(r > r= 



E 



i=q,q,g 



X dx 



(A.12) 



At NLO the thrust cut removes any regions of phase-space involving divergences, so we may 
set e = 0. Only the real radiation contributes, yielding 



ddjj 
dz 



dog 
dz 



(r >T'') = a 



(0) 



a 



27r 

(0) OisCp 



(r > r'^) = a 



(0) 



asCp 



a 



27r 

(0) OisCp 
TT 



e{l - z - T^)9{z - 2t^) 
9{1 - z - T'')eiz - 2t'') 

d{l- Z-T^)e{z-2T'') 



dx 



+ x^ 



-z+T- {'^-Z){l-X) 



rl„ , z-t" z2 - 2z(r= + 2) + Sr'^ 
-Pgg(z)ln — + 



dx 



4(1 - z) 

(2 - X - z)2 + 



l-z+T- \X + Z-I){l-X) 

C 

e{l-z- T^)e{z - 2r^) Pgq{z) In -z + 2t^ 



(A.13) 



A. 4 NLO singular vs. nonsingular contribution 

The singular cross section can directly be obtained from our resummed calculation in eq. (2.2), 
by choosing the renormalization scale of all the objects equal to /x, which switches off the 
resummation. By subtracting this from the NLO cross section calculated in the previous 
section, we obtain the nonsingular piece. In fig. 8 the NLO cross section for z = 0.2 and 
z = 0.8 is shown separated into its singular and nonsingular contribution. A first cross-check 
of our calculation is provided by the fact that the nonsingular piece is finite in the r — )• 
limit, whereas both the full NLO and the singular piece diverge. The contribution of the 
nonsingular to the total cross section is suppressed for small r by 0{t). However, for larger 
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values of r its contribution becomes more important, and above the NLO kinematic threshold 
r > min{l/3, ^ — z} the total cross section vanishes and the singular and nonsingular exactly 
cancel each other. It is in this region that the inclusion of the nonsingular is essential to 
avoid a negative differential cross section in r or an unphysical turn-over in the cross section 
integrated up to r < r'^. 



B Choice of running scales 



The resummation of logarithms of r and 1 — z, as well as large vr^ terms, is accomplished 
by evaluating H, J, J and S in eq. (2.2) at their natural scales /xj, /ij-, and using 
their respective RGEs to evolve them to a common scale /U. Having made some improvements 
compared to our earlier work ref. [20], we discuss this in some detail here. 

The hard function K is the square of a time-like form factor and contains large vr^-terms 
at [III = Q arising from \v?{—\Q/pLii)- These can be resummed through the complex scale 
setting nil = —'\Q [41-44]. Following ref. [32], we note that there are three distinct kinematic 
regions where the resummation of the logarithms of r must be handled differently: 



1) m ^ -iQ , fJ-J - \/AqcdQ , 

2) fiH =i -iQ , fJ-J - VrQ , 

3) ifiH = fJ-j = fJ-s - Q ■ 



= Aqcd 



Region 1) is not of interest for our application, and requires a completely nonperturbative 
soft function. The resummation of the threshold logarithms of 1 — 2: sets the scale fi j different 
from /xj for the diagonal terms Jqq [22]. In addition, the measurement of z affects the upper 
bound T < Tmax = min{0.5, 1 — z}, implying that region 3) sets in at smaller values of r for 
z > 0.5. Taking this into account, a smooth transition between the three regions is given by 



r 1 



Hj{t, z)= l + ej 6'(r3rmax - 

l^ji'T, = Vl- Z fij{T, Z) , 

IJ'S{t,z)= 1 + 65 6'(r3rmax - t) ( 1 



T 



\/\P'H\ /^run(T"/Tmax, Z, 



/^mnlT/Tmax, Z, \^Xh\) 



(B.l) 



where scale uncertainties will be estimated by varying en, ej and as shown in table 1. In 
addition to the explicit 2;-dependence in fij, there is an implicit dependence through Tj^^x* 
For the profile ^J-run{T/Tma.x, z, l/ii^l) we use a combination of two quadratic functions and a 
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var 


gh 






interpretation 


1 


-1 
i 








central value 


2 


2 














hard scale variation 


3 


0.5 










4 


1 


0.5 





jet scale variation 




1 


5 


-0.5 





6 


1 





0.5 


soft scale variation 




7 


1 





-0.5 





Table 1. The scale variations in terms of the parameters in eq. (B.f ). 
linear function as in refs. [32, 34], 



Ho + ar^/ri r < ri , 

2ar + b ri < r < r2 ■ 

\iJ-h\ - a{r - r3)'^/{r3 - r2) r2<r<r^. 

, \^^H\n - ^^o{r2 + n) 



(B.2) 



ri-r2- rs ri - r2 - 

The expressions for a and h follow from demanding that //run is continuous and has a con- 
tinuous derivative. The value of /io determines the scales at r = 0, while ri^2,3 determine 
the transition between the regions discussed above. For r > raTmax, our choice for //run en- 
sures that the resummation of logarithms of r turns off (but vr^ resummation and threshold 
resummation are still present). 

For the cross section integrated up to r'^, we use the same scale choices as in eq. (B.l) 
with r — )• r"^. In moment space, Tmax = 0.5 and the equation for fij is modified to 



(B.3) 



We complete our description by specifying our choice for the parameters of the profile function 
in eq. (B.2). Our central value corresponds to: 

4GeV 



en = 1 , ej = 65 = . 



2GeV. 



ri 



Q 



r2 = 0.5 , rs = 1 . 



(B.4) 



As mentioned before, we estimate the perturbative uncertainties by a combination of separate 
scale variations, which are listed in table 1. In our cross section plots, we show the envelope 
of these variations. In the tables in app. C we give results for all the separate scale variations, 
allowing one to take correlations in the uncertainties into account. 

C Perturbative coefficients Cj in moment space 

In this appendix we give numerical results for the coefficients in eq. (4.1) at NNLL+NLO. 
Tables 2 and 3 contain the results for Cq and Cg, respectively. For each of the independent 
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scale variations in table 1 there is a separate set of Cj, labelled by the number in the first 
column. By extracting the Dj{N,fi = 1 GeV) for each of these variations, we obtain an 
estimate of the perturbative uncertainties for the FFs, including correlations between these 
uncertainties. 
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1.65 
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1.78 


1.77 
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1.95 


1.89 


1.85 


1.82 


1.80 


1.79 


5 


0.15 


1.83 


1.85 


1.86 


1.88 


1.88 


1.89 




0.16 


1.89 


1.89 


1.89 


1.89 


1.88 


1.88 




0.17 


1.94 


1.92 


1.90 


1.89 


1.88 


1.87 
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1.98 


1.94 
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1.88 
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1.88 
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2.05 
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1.88 


1.85 


6 


0.15 


1.84 


1.85 


1.85 


1.86 


1.87 


1.88 




0.16 


1.91 


1.89 


1.89 


1.88 


1.88 


1.89 




0.17 


1.96 


1.93 


1.91 


1.90 


1.89 


1.88 




0.18 


2.00 


1.96 


1.93 


1.91 


1.89 


1.88 




0.19 


2.04 


1.98 


1.94 


1.92 


1.89 


1.88 




0.20 


2.07 


2.00 


1.96 


1.92 


1.89 


1.87 


7 


0.15 


1.48 


1.50 


1.52 


1.54 


1.57 


1.59 




0.16 


1.59 


1.60 


1.60 


1.62 


1.63 


1.65 




0.17 


1.69 


1.68 


1.67 


1.68 


1.68 


1.69 




0.18 


1.77 


1.74 


1.73 


1.72 


1.72 


1.72 




0.19 


1.83 


1.80 


1.77 


1.76 


1.75 


1.75 




0.20 


1.89 


1.84 


1.81 


1.79 


1.77 


1.76 



Table 2. Cg(r^ N,Q = 10.52 GeV, /i = 1 GeV) in units of nb at NNLL+NLO, with scale variations. 
Cq applies to both up and down-type quarks and was defined in eq. (4.3). 
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057 
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1 85 
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0.18 
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0.382 
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0.409 


0.246 
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0.094 
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0.486 


0.307 


0.212 
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0.099 




0.17 


0.526 


0.328 


0.225 


0.166 


0.129 


0.104 




0.18 


0.564 


0.346 


0.236 


0.174 


0.134 


0.108 




0.19 


0.598 


0.364 


0.246 


0.180 


0.139 


0.111 




0.20 


0.630 


0.379 


0.255 


0.186 


0.143 


0.115 


O 


1 5 




01 9 


-0 01 1 


-0 030 


-0 043 


-0 053 




Ifi 


f) 083 
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-0 004 
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-0 039 


-0 049 
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0.042 
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-0.019 


-0.034 


-0.044 
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0.128 


0.054 


0.011 


-0.013 


-0.029 


-0.040 




0.19 


0.149 


0.065 
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-0.007 


-0.024 


-0.035 
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-0.030 


4 


1 ,5 
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1 fi 
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095 


075 
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0.431 


0.265 
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0.130 


0.099 


0.078 
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0.457 


0.278 


0.186 


0.135 


0.103 


0.081 




0.19 


0.481 


0.289 


0.193 


0.139 


0.106 


0.083 




0.20 


0.503 


0.300 


0.200 


0.144 


0.109 


0.085 


(J 


1 5 


n 086 


045 


n 01 2 


-0 01 


-0 025 


-0 036 




0.16 


0.131 


0.073 


0.032 


0.007 


-0.010 


-0.022 




0.17 


0.174 


0.099 


0.051 


0.023 


0.004 


-0.008 




0.18 


0.215 


0.123 


0.069 


0.037 


0.017 


0.003 




0.19 


0.254 


0.146 


0.086 


0.051 


0.029 


0.014 




0.20 


0.291 


0.168 


0.101 


0.063 


0.039 


0.023 


6 


0.15 


0.272 


0.178 


0.120 


0.084 


0.060 


0.044 




0.16 


0.306 


0.195 


0.130 


0.091 


0.066 


0.048 




0.17 


0.338 


0.211 


0.140 


0.098 


0.071 


0.053 




0.18 


0.368 


0.227 


0.149 


0.104 


0.076 


0.057 




0.19 


0.396 


0.241 


0.158 


0.111 


0.081 


0.061 




0.20 


0.422 


0.254 


0.166 


0.117 


0.086 


0.065 
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1 5 


n 236 


1 53 


1 03 
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053 
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16 


272 


1 72 


115 

V7. -L -L tJ 


081 


059 
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\) . \J tit: 




0.17 


0.305 


0.190 


0.127 


0.089 


0.066 


0.049 




0.18 


0.337 


0.207 


0.137 


0.097 


0.071 


0.054 




0.19 


0.368 


0.223 


0.147 


0.104 


0.077 


0.058 




0.20 


0.396 


0.239 


0.157 


0.111 


0.082 


0.062 



Table 3. (^^(t^ N,Q = 10.52 GeV, = 1 GeV) in units of nb at NNLL-hNLO, with scale variations. 
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